Evidence of Orbital Motion in the Binary Brown Dwarf Kelu-lAB 
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ABSTRACT 

We have resolved Kelu-1 into a binary system with a separation of ~290 mas 
using the Laser Guide Star Adaptive Optics system on the Keck II telescope. 
We have also re-analyzed a 1998 HST observation of Kelu-1 and find that the 
observed PSF is best fit by a binary object separated by 45 mas. Observations 
on multiple epochs confirm the two objects share a common proper motion and 
clearly demonstrate the first evidence of orbital motion. Kelu-IB is fainter than 
Kelu-lA by 0.39±0.01 magnitudes in the K' filter and 0.50±0.01 magnitudes in 
the H filter. We derive spectral types of L2±l and L3.5±l for Kelu-lA and B, 
respectively. The separation of flux into the two components rectifies Kelu-1 's 
over-luminosity problem that has been known for quite some time. Given the 
available data we are able to constrain the inclination of the system to >81° and 
the orbital period to >40 years. 

Subject headings: stars: individual (Kelu-lAB) — stars: low-mass, brown dwarfs 
— techniques: high angular resolution 



1. Introduction 



L dwarfs are a class of objects with spectral types cooler than the latest M dwarfs 
(Kirkpatrick et al. 1999). Their population consists of both substellar (i.e. brown dwarfs) 
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and stellar mass objects, with the brown dwarfs being identified by the presence of Li I 
absorption at 6708 Ain their spectra. Kehi-1 was one of the first sTich brown dwarfs discovered. 
It was found at a trigonometric distance of 18.7 pc (Dahn ct al. 2002) as part of a proper motion 
survey (Ruiz, Leggett, k. Allard 1997). Since Kelu-l's discovery, more than 400 L dwarfs have 
been classified -"^ in large surveys such as the DEep Near-Infrared Survey (DENIS), the 2-Micron 
All-Sky Survey (2MASS), and the Sloan Digital Sky Survey (SDSS). A spectral class sequence 
was established (Kirkpatrick et al. 1999, 2000; Martin et al. 1999b), and Kelu-1 was subsequently 
assigned a spectral type of L2 in the optical (Kirkpatrick et al. 1999) and L3 in the near-IR (Knapp 
et al. 2004). These other discoveries provided a baseline against which to compare Kelu-l's other 
characteristics. For example, it was quickly realized that Kelu-1 was much brighter than other L2 
dwarfs, causing Martin et al. (1999b) to first suggest Kelu-1 a binary. This super-luminosity has 
been seen in other studies (Leggett et al. 2001; Golimowski et al. 2004) where the main source of 
this over-brightness was attributed to either an unresolved binary or a young age. 

Martin, Brandner, & Basri (1999a) observed Kelu-1 on 1998 August 14 with the NICMOS 
camera on HST and no companions were found for separations greater than 300 mas and magnitude 
differences <6.7 mag (F165M filter). Koerner et al. (1999) observed Kelu-1 in a non-AO, seeing- 
limited program on Keck I using the NIRC instrument on 1998 February 14 and 1999 February 
9. Because Kelu-1 had a fairly circular point spread function during these observations (PSF 
FWHM!=a0.27"), it was used as a PSF star for the other companion search targets of their project. 
While these observations did not discount Kelu-1 as being a very close, unresolved binary, they did 
warrant a look at the young age hypothesis. 

Based on the strength of the Li I absorption feature, Basri et al. (1998) constrained the 
age of Kelu-I to be 0.3-1 Gyr. With this age range and Kelu-l's brightness, it should have an 
effective temperature of 2100-2350 K, about 400 K higher than objects with similar spectral types 
(Golimowski et al. 2004). For it to have an effective temperature consistent with other L3 dwarfs, 
Golimowski et al. indicate Kelu-1 must have an age of -^10 Myr, more than an order of magnitude 
lower than the Li I-established lower age limit. These discrepancies mean that a young age is 
unlikely to be the cause of the over-luminosity. 

Aside from the luminosity, there is other evidence that supports the binarity of Kelu-1. It 
is reported to be photometrically and spectroscopically variable with a period of 1.8 hrs (Clarke, 
Tinney, k, Govey 2002; Clarke, Tinney, & Hodgkin 2003). If this period corresponds to an orbital 
period, then the system would likely be circularized and the orbital radius of the system would be 
less than a solar radius (Clarke et al. 2002). In addition, Clarke et al. (2002) claim that the ~120 
km s~^ radial velocity produced by a very close binary such as this could account for the large line 
broadening (60ib5 km s~^) measured by Basri et al. (2000). 

We are conducting an adaptive optics (AO) survey of L and T dwarfs at Keck Observatory in 



See iittp://dwarfarchives.org for a complete list of all known spectroscopically determined L dwarfs. 
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an effort to determine binarity in these objects (an article describing the full details of this project 
is in preparation). Because of the intrinsic faintness of L and T dwarfs (magnitudes for Kelu-1 arc 
V=22.1 and R=19.2; Ruiz ct al. 1997), it is not possible to use them as Natural Guide Stars (NGS) 
where the magnitude limit is about y=14. Instead, we utilize the newly commissioned Laser Guide 
Star Adaptive Optics system (LGS AO)^ where the wavefront corrections are performed on a layer 
of sodium atoms 80 km high in the Earth's mesosphere that are excited by a 10-14 Watt laser 
tuned to 598nm. The projected laser and returned light propagate along the same path, requiring 
the use of a natural star (or other point-like astronomical source) for tip-tilt measurements. The 
tip-tilt reference can be as faint as i?=18.5, but must be within 60" of the science target. Thus, the 
area of sky reachable with AO is greatly increased by the less stringent requirements for a reference 
star provided by the used of the LGS, and many L and T dwarfs binaries can be observed with 
high resolution imaging from the ground. We included Kelu-1 in our survey in an effort to resolve 
the controversy over its binarity and over-luminosity and to search for any low-mass companions. 
In this article we discuss the detection of a companion to Kelu-1, first reported in Liu &i Leggett 
(2005), and show the first evidence of orbital motion of the companion. 



2. Observations and Data Reduction 

2.1. Keck LGS AO 

Observations of Kelu-1 were obtained with NIRC2 (K. Matthews et al. in preparation) behind 
the LGS AO system on 2005 March 4 and 2005 April 30 during shared-risk time. A suitable tip- 
tilt reference star (ID=0643-0289866, i?=14.3, separation=27") was identified in the USNO-B 1.0 

catalog (Monet et al. 2003). Skies were photometric on both nights and the AO system provided 
very good corrections on our targets. Kelu-1 was observed with the K' filter during both epochs, 
with additional exposures in the H filter during follow-up observations in April. These filters were 
chosen to maximize the Strehl ratio in our images. The data were reduced using custom and public 
IDL scripts. 

For each night, two sets of reduced images were produced from the data: a photometric set 
and an astrometric set. The only difference in these sets is that an image distortion correction was 
performed on the astrometric data set. The distortion in the narrow camera is negligible within 
a region a few hundred pixels on a side around the center of the chip and gradually increases to 
as much as 4-5 pixels near the edge. While the correction fixes the distortion with residuals of 
0.81 and 0.62 pixels in the x and y directions, respectively, it does not preserve flux information. 
Consequently, its use in images for photometric analysis is not suitable. 

The plate scale for both nights was determined using astrometric images of the core of the 



^http: / / www2.keck.hawaii.edu/ optics/lgsao 
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globular cluster M5. This field was observed with WFPC2 on HST, with the stars in common with 
our observations present in the Planetary Camera. We imaged this field with a 4 (in April) or 
5 (in March) position dither pattern in the narrow camera (~10" field-of-view). From the WCS 
information in the HST image headers, we determined the sky positions of several stars present in 
our NIRC2 images. These positions were then used to calculate the plate scale and rotation angle 
for each position in the dither. The average of these values is the resultant scale and rotation for 
that night. The errors in these quantities are equal to the standard deviation of the measurements. 
Table 1 presents the plate scale for these nights. 

In our March 4 images of the Kelu-1 field, we saw what appeared to be a binary object with 
a separation '--'300 mas and a differential K' magnitude of ~0.4. In order to obtain proper mo- 
tion confirmation of the companion, we re-observed this system on April 30. On both epochs, we 
determined the separation and position angle of the candidate companion by using an average of 
measurements on the individual reduced astrometric images. The errors are the standard devia- 
tions of these measurements. Table 2 shows these results. The moderately large proper motion of 
Kelu-1 (285.0±1.0 mas yr"^ at 272.2±0.2°; Dahn et al. 2002) made it possible to confidently con- 
firm physical companionship based on observations separated by only two months. The measured 
positions of the companion were about 15a away from its predicted background position, clearly 
demonstrating that these two objects have a common proper motion and are physically associated. 
Our binary conclusion firmly supports the recent announcement of Kelu-l's binarity from Keck 
LGS AO observations obtained in May 2005 (Liu k Leggett 2005). 

For both nights, all of the photometric images were offset such that northeast object was in the 
same location in the overlap region. This image stack was then median averaged to obtain the final 
image (Figure 1), on which the photometry was performed. The components were well separated 
and no PSF subtraction routines were needed. No photometric standards were observed on either 
night, so only differential photometry is shown in Tabic 3. The northeast object is the brighter 
component in both the H and K'. It will henceforth be called Kelu-IA and the fainter southwest 
component Kelu-IB. Our measured astrometry and differential photometry are consistent with 
those of Liu & Leggett (2005). 

2.2. HST NICMOS 

We re-examined E. Martin's 1998 HST NICMOS images using the PSF modeling technique 
described in Stephens &; Noll (2006). In short, this technique involves the sub-pixel shifting of two 
model PSFs and iteratively changing the relative brightness of the PSFs to obtain the minimum 

residual for each shifted pair. The lowest minimum residual for all shifted pairs represents the 
best fit for the observed PSF. Our analysis discovered that the observed PSFs in all three filters 
(FllOM, F145M, and F165M) were better fit by double PSF solutions than they were by single 
PSF solutions. The derived astrometry of the components is presented in Table 2. Because the 
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objects are unresolved, a 180° ambiguity exists in the value of the position angle. 

It is also possible to estimate the magnitudes of the two components (listed in Table 4) using 
this technique. Because the two components are unresolved, the determination of the individual 
magnitudes is subject to systematic uncertainties in the fitting process. Therefore, it is difficult to 
quantify the uncertainty in these magnitudes. Consequently, although these calculations may not 
determine the actual difference in magnitudes between the two components, they do indicate the 
reliability of a binary solution. Furthermore, the computed magnitudes are also correlated, so if one 
object is actually fainter than computed, the other component would be brighter, thus conserving 
the total flux in the system. 

3. Discussion 
3.1. Component Characteristics 

The spectral types of the components can be estimated using the photometry and absolute 
magnitude-spectral type relation in Vrba et al. (2004) and the differential photometry presented 
here. For this exercise we assume that the AiJ and AX' measurements in this study are equivalent 
to Ai? and Aif in the CIT photometric system as employed by Vrba et al. While we know that 
this is not the case, our interest is only in approximately determining the spectral types of the 
components; the only way to determine spectral types for these or any other objects is to obtain 
properly calibrated spectra. Based on the multi-system, near-IR photometry study of Stephens & 
Leggett (2004), wc conservatively impose an additional photometric error of 0.15 mag for using this 
assumption. This error is added in quadrature to our measured photometric errors. 

The i^ciT and i^ciT apparent magnitudes of Kelu-lAB are 12.46±0.04 and 11.80±0.03, respec- 
tively (Dahn et al. 2002). Coupled with our measured Ai7=0.50±0.15 and A/s:'=0.39±0.15, the 
apparent magnitudes of the resolved components become i?[A]=12.99, K[A]=12.37, i7[B]=13.49, 
and if[B]=12.76 (Table 4). Converting these magnitudes from apparent to absolute and applying 
equations (2) and (3) from Vrba et al. (2004) yield a spectral type of ~L2±1 for Kelu-IA and 
~L3.5ibl for Kclu-IB. These estimates are consistent with the optical and near-IR spectral types 
of the composite system and with the spectral types derived by Liu &; Leggett (2005). 

With this photometry and spectral type estimates, we can place these objects on magnitude- 
spectral type diagram. Figure 3 presents a modified version of Figure 4 from Vrba et al. (2004). In 
agreement with Liu & Leggett (2005), it is quite clear that the over-luminosity of Kelu-lAB was 
caused by the unresolved binary and not due to an unusually young effective temperature-based 
age (Golimowski et al. 2004). 

Finally, we can estimate the effective temperature from the photometry presented here. Golimowski 
et al. (2004) calculate a if -band bolometric correction (BCi^-) for Kelu-lAB of 3.31±0.07. While 
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this calculation is independent of the spectral type of Kelu-lAB, the spectral type used for Kelu- 
lAB and all of the other objects in that study were based on near-IR spectra. The spectral types 
for Kelu-IA and B presented here are based on a relation derived from optical spectral types. The 
near-IR spectral type of Kelu-lAB is L3±l (Knapp et al. 2004), whereas the optical spectral type 
is L2±0.5 (Kirkpatrick et al. 1999). As discussed in Kirkpatrick (2005), it is neither uncommon nor 
unexpected for L and T dwarfs to be classified differently in the optical and near-IR. Consequently, 
it is unclear what value of BCk should be used from Golimowski et al. (2004). 

Fortunately, for early L dwarfs such as Kelu-1 A and B, the range of BC^^ is quite small. 
The average value of BC^^ for objects with near-IR spectral types L0-L5, inclusive, is 3.32 with a 

standard deviation of 0.07. Since this average value is consistent with the value of BCx computed 
for the unresolved Kelu-lAB, we choose the unresolved value as the bolometric correction for 
the individual components. We note, however, that the uncertainty in the effective temperature 
determination is largely dominated by the unknown age of Kelu-1 and, given the constraints on the 
bolometric corrections of early L dwarfs, the exact value of BC^ is unimportant. 

Based on our photometry, the absolute K magnitudes of Kelu-IA and B are 11.02 and 11.41, 
respectively. Assuming a if -band bolometric correction of 3.31 for each component yields bolomet- 
ric magnitudes of 14.33 and 14.72 and luminosities (logio L/Lq) of -3.83 and -3.99. Finally, using 
an age range of 0.3-1 Gyr (Basri et al. 1998) and the Chabrier et al. (2000) DUSTY atmosphere 
models, we determine an effective temperature range for Kelu-IA of 1900-2100 K and for Kelu-IB 
1700-1900 K (Table 4). Liu & Leggett (2005) computed effective temperatures for Kelu-IA and B 
(2020 K and 1840 K) by scaling the effective temperature of Kelu-lAB (Gohmowski et al. 2004) by 
their luminosity calculations. Despite the differences between our method and theirs, our results 
are completely consistent. 

3.2. System Chciracteristics 

It is not known if the Li I absorption at 6708 A detected by Ruiz et al. (1997) and others is 
from one or both components. Liu & Leggett (2005) suggest that the presence of any lithium in 
the unresolved spectrum indicates both components are substellar and, at the very least, Kelu-IB 
bears lithium in its atmosphere. Consequently, the masses of both components must be <0.065 
Mq (Chabrier &: Baraffe 1997; Basri 1998). To estimate component masses, we use the photometry 
derived here along with the age estimate of Basri et al. (1998) to place the system components 
on the DUSTY atmosphere model evolutionary tracks (Chabrier et al. 2000). Prom the combined 
mass estimates presented in Figure 2, we estimate the mass of Kelu-IA to be 0.060±0.01 Mq and 
that of Kelu-IB to be 0.055±0.01 M©, in good agreement with the mass estimates derived by Liu 
& Leggett (2005). Even though the conservative errors on these masses do include stellar masses, 
as mentioned above the true masses are likely to not be greater than about 0.065 Mq. Of course, 
by following the orbital motion of the companion over time, it will be possible to derive dynamical 
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masses, which will provide much needed constraints on the models. 

The computation of these dynamical masses will be some time in coming. Table 2 presents 
all of the known astrometry of the Kelu-lAB system as of 2006 January. It is quite clear from this 
astrometry that Kelu-IB has exhibited orbital motion from the time it was first resolved (our 2005 
March 4 observation) to the time of its last observation (a public archival HST observation on 2005 
July 31 obtained by W. Brandner). In fact, the displayed motion suggests that the companion is 
still moving toward apoapsis. During this nearly 5 month period the rate of change in the separation 
has not changed significantly and there is insufficient phase coverage to determine a reliable orbit. 
Nonetheless, we can make some crude estimations of the orbit based on the available data. 

Given so few points in the orbit, it is impossible to determine how much of the nearly linear 
motion is caused by a high inclination, very high eccentricity, or both. In order to simplify the 
discussion, we assume that the orbit is circular and thus, is being viewed nearly edge-on. 

The most recently measured 2005 positions and the 1998 PSF-modeled position provide some 
tight constraints on the inclination of the orbit. Although the position angle of the 4 resolved 

measurements is not significantly changing, it is not possible to extend a line connecting all four of 
these points directly back to Kclu-IA. Therefore, the inclination of the orbit must be <90°. We can 
get a rough estimate of the minimum inclination by assuming the maximum extent perpendicular to 
the 2005 July position is the separation of the modeled PSFs in 1998, i.e. 45 mas. This assumption 
results in an inclination >81° and is much more constrained than the minimum inclination suggested 
by Liu &; Leggett (2005) based on the size of the NICMOS PSF. Finally, since it is quite clear that 
Kclu-IB is still moving away from Kelu-IA, this minimum inclination estimate will only increase 
with time. 

The period of the binary can be estimated using the 2005 July 31 separation (298 mas or 5.4 
AU at the distance of Kelu-1) and the masses derived above. With the caveat that the most recent 
separation is not the maximum orbital separation, the calculated period is > 37 years. This period 
is important for two reasons. First, it means that it is possible to obtain dynamical masses for the 
components in a reasonable time-frame, something that has been done for only a few other binary 
brown dwarf systems (Brandner et al. 2004; Zapatero Osorio et al. 2004; Bouy et al. 2004). Second, 
this period means that the 1.8 hour period (or 3.6 hours if the variations are ellipsoidal in nature) 
reported by Clarke et al. (2002, 2003) is not due to this binary. Perhaps one of the components is 
itself a very close binary, or perhaps the shorter period is simply the rotation period for one of the 
components. In addition, the 60 km s~^ rotation velocity measured by Basri et al. (2000) on 1997 
June 2 is from a composite spectrum. The maximum orbital velocity for Kelu-lAB is only ~3 km 
s~^, which translates to a displacement of ib0.09A and ibO.OSA for the Cs I and Rb I atomics 
lines used by Basri et al. This wavelength shift cannot account for the large equivalent widths of 
these lines (1.7A and 2.54A for Cs I and Rb I, respectively), so some other mechanism must be at 
work. In time it might be possible to separate the lines from each component using high resolution 
spectroscopy and resolve some of these issues. 
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4. Conclusions 

We have used the LGS AO system on the Keck II telescope to observe the binary brown dwarf 
Kelu-IAB. We have also re-examined the 1998 HST observations of Kelu-1 and show that the PSFs 
in those images are best fit by a binary object solution. Images from multiple epochs confirm the 
pair as having a common proper motion and demonstrate that Kelu-IB is still heading towards 
apoapsis. Wc constrain the inclination of a circular orbit to >81°. While this companion detection 
does rectify the over-luminosity "problem," it cannot account for the 1.8 hour photometric period 
or the 60 km s~^ rotation velocity. Periodic monitoring of this system with high resolution imaging 
and spectroscopy is needed to ensure the maximum extent of the orbit is observed and measured, 
allowing for a more robust determination of the physical orbit. 
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operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract 
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University of California and the National Aeronautics and Space Administration. The Observatory 
was made possible by the generous financial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant cultural role and reverence that the summit 
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Table 1. NIRC2 Plate Scales and Rotation Angles 



Date 


Scale 


Rotation 


Dither 


(UT) 


(mas/pixel) 


(deg) 


Positions 


2005 Mar 4 


9.97±0.02 


+0.08±0.04 


5 


2005 Apr 30 


9.96±0.02 


+0.11±0.03 


4 



Table 2. Kelu-IB Astrometry 



Date 


Separation 


RA. 


Reference 


(UT) 


(mas) 


(deg) 




1998 Aug 14 


45±18 [ 


38.0 or 218.0]±11.9^ 


this work 


2005 Mar 4 


284.3±0.8 


220.9±0.3 


this work 


2005 Apr 30 


289.3±0.7 


221.1±0.1 


this work 


2005 May 1 


291±2 


221.2±0.6 


Liu & Leggett (2005) 


2005 Jul 31 


298±3 


221.3±0.9 


this work^ 



'^Because the objects are unresolved, there is a 180° ambiguity in the 
position angle solution. 

^Based on HST observations by W. Brandner that are present in the 
public archive 



Table 3. Kelu-lAB Photometry 



Date (UT) Filter Amag 



2005 Mar 4 
2005 Apr 30 



K' 
K' 
H 



0.39±0.01 
0.39±0.01 
0.50±0.01 
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Table 4. Kelu-lAB Derived Properties 



Property 


A 


B 


Spectral Type 


L2±l 


L3.5±l 


H magnitude 


12.99±0.15 


13.49±0.15 


K magnitude 


12.37±0.16 


12.76±0.16 


FllOM magnitude 


14.30±0.04 


15.13±0.10 


F145M magnitude 


13.83±0.10 


14.04±0.13 


F165M magnitude 


13.03±0.07 


13.35±0.11 


Mass (Mq) 


0.060±0.01 


0.055±0.01 


Te// (K) 


1900-2100 


1700-1900 
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Fig. 1. — K' image of Kelu-lAB taken on 2005 March 4. The image is 1.2"on a side with 
North up and East to the left. Kelu-IB is located 284.3±0.8 mas at a position angle of 
220.9±0.3 deg from Kelu-IA. 
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Fig. 2. — a. Absolute H magnitude as a function of age for Kelu-lA, Kelu-IB (boxed regions) 
and the DUSTY models (dotted lines) for masses from 0.1 (top) to 0.02 M©. b. Same 
as a except for absolute K magnitude. The age range is derived from the strength of the Li 
1 absorption feature in the composite spectrum Basri et al. (1998). The combined estimate 
for the masses of Kelu-1 A and B are 0.060±0.01 and 0.055±0.01 M©, respectively. 
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Fig. 3. — Absolute K magnitude as a function of spectral type for the objects in Table 
7 of Vrba et al. (2004). The position of the unresolved Kelu-lAB is marked as an open 
square; the individual components are shown with solid squares. The binarity rectifies the 
well-documented over-luminosity in the K band. 
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